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1 INTRODUCTION 



The halo of our galaxy is believed to have formed from the 
gravitational collapse of dark matter (DM) and a continuous 
merger with other haloes. The formation of haloes, described 
by the Jeans- Vlasov-Poisson equation, proceeds with the for- 
mation o f dark matter sh ells, a lso referred to as d ark matter 
caustics JSikivie fc Ipser 1992; Sikivie et al. Ill997l:[ 
Colombill2005i : iNataraia 



ABSTRACT 

Gravitational collapse of dark matter, merger of dark matter haloes and tidal disrup- 
tion of satellites are among processes which lead to the formation of fine and dense 
dark matter shells, also known as dark matter caustics. The putative weakly inter- 
acting species which may form the dark matter are expected to strongly annihilate in 
these dense regions of the Milky Way halo and generate in particular antiprotons and 
positrons. We derive the flux of these rare antimatter particles at the Earth and show 
that it depends significantly on the cut-off radius of the dark matter distribution at the 
galactic centre. Boost factors of 30 are found with respect to a smooth NFW profile 
for high-energy antiprotons and low-energy positrons if this cut-off radius is taken to 
be 300 pc - a somewhat extreme value though. This yields a detectable antiproton 
signal around hundreds of Gev in models where the annihilation cross section today 
is enhanced by non-perturbative effects as in the generic case of a heavy Wino. How- 
ever, dark matter caustics cannot provide a better explanation for the HEAT excess 
reported above 10 GcV than a smooth NFW or isothermal cored distribution. 

preprint no: LAPTH-1235/08 

rounded by a large number of primordial shells and fossil 
shells of disrupted satellites. The density of the shells can 
exceed that of the diffuse background by an amount which 
depends on the age of the galaxy and the number of times 
the satellite has wrapped around our galaxy. Due to their 
vicinity and over-density, these shells can be important for 
dark matter search experiments. 
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the density peaks collapse under self-gravity, at the surfaces 
of (formally) infinite density where dark matter streams 
meet, caustics emerge. The infalling satellites are also dis- 
trupted by a mechanism similar to the formation of the pri- 
mordial caustics. As a satellite moves inside the potential 
well of its host halo, a tidal tail forms, and around the apap- 
sis of the orbit high-d ensity shells or caustics emerge [see e.g. 
iHavashi et al] (|2003l )]. Since each infalling satellite also has 
his own caustics, the hierarchical formation of structures 
yields a hierarchy of caustics. Observational examples for 
the caustics are the shell galaxies, where shells of stars form 
during the merger of galaxies, by a mechan ism very similar 
to the formation of dark matter caustics (iMalin fc Carted 
1 19801 : ICarter et al.]|l982l : iHernguist fc Quinnll 19881 , 1 19891 ) . 

Once formed, they cannot be destroyed : caustics are 
permanent structures in real and in phase space. However, 
their over-density with respect to the background density 
can diminish with time due to continuous increase in the 
number of streams. Therefore, it is likely that we are sur- 



The nature of dark matter remains a mystery. Super- 
symmetry and extra-dimension extensions of the standard 
electroweak model provide a natural candidate in the form of 
a weakly interacting and massive particle (hereafter WIMP). 
These species should fill up the galactic halo. If dark mat- 
ter consists of WIMPS then they are expected to strongly 
annihilate in the dense regions of our halo (i.e. in the caus- 
tics) and generate in particular gamma-rays and charged 
cosmic rays such as positrons and antiprotons. Indeed, due 
to their very short diffusion length, positrons can be excel- 
lent tracers of nearby caustics. Although the implication of 
substructures, or surviving satellites, in boosti ng the cosmic- 
ray signal has been st udied extensively [see e.q. lLavalle et al] 
(|2007al ): lLavalle et al., (2007bi) ]. the signature of dark matter 
caustics on cosmic-rays has been rarely considered. 

Here, we propose to study the boost in positron and an- 
tiproton signals due to annihilation in dark matter caustics. 
The ideal way to pursue this study would be to use a very 
high-resolution simulation. However, even the latest simu- 
lations have not yet achieved a high enough resolution to 
identify the caustics in a ty pical galaxy, although imp ressive 
progresses are being made (jVogelsberger et al.l 120071 ). Ana- 
lytic models for the formation of caustics are mainly based 
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on the self-similar seco ndary infall (|Fillmore &i GoldreichI 
1 19841 : iBertschingedl 1985h . We use a generalized version of the 
secondary infall model whic h takes into account the finit e 
velocity dispersion of DM (jMohavaee fc ShandarinI l2006h . 
Although secondary infall model has serious limitations be- 
cause of the assumption of spherical symmetry and smooth 
accretion and because it ignores the hierarchical scenario of 
structure formation, it has so far provided a paradigm for 
the st udy of dark matter haloes [e.g. see lAscasibar et alJ 
l|2007l )]. 

In Section 2, we review and generalize the secondary 
infall model for dark matter with finite velocity dispersion 
and give an analytic fit for the density profile which includes 
the caustics for a typical galaxy like our own. The cosmic 
ray signature of dark matter caustics is studied in Section 3 
for large antiproton and positron horizon radii. In Section 4, 
the possibility of detecting a nearby caustic through its short 
range positron signal is investigated. We finally conclude in 
Section 5. 



2 THE SECONDARY INFALL OR 

SELF-SIMILAR ACCRETION MODEL 

Analytic evaluation of the halo density profile, and the pre- 
diction of the existence of caust i cs insi de th ese str u ctures , 
started with the works of iGottI l\197^ and ICunnl (|l977t ). 
With the main objective of explaining the flattening of the 
rotation curves of the galaxies, they considered the forma- 
tion of a dark matter halo from the infall of matter onto 
an already-formed galaxy (or in later works onto a spher- 
ical overdense region). In an Einstein-de Sitter Universe a 
spherical over-density expands and then turns around to col- 
lapse. After collapse and at late times, the fluid motion be- 
comes self-similar : its form remains unchanged when length 
are re-scaled in terms of vta, the radius of the shell which 
is currently turning around and collapsing. Physically self- 
similarity arises because gravity is scale-free and because 
mass shells outside the initial over-density are also bound 
and turnaround at successively later times. Self-similar so- 
lutions give power-law density profiles whose exact scal- 
ing properties depend on the initial condition, the central 
boundary conditi ons and on whether the fluid is collision- 
less o r coUisional (|Fillmore fc GoldreichI [T984I : iBertschingeij 
1 19851 ). The global density profile is a power-law which for 
suitable choices of initial conditions provides an explanation 
of the flattening of the rotation curves of the galaxies. A den- 
sity proflle of results for the galactic-scale fluctuations. 
However, the density profile is convolved with many spikes 
{i.e. caustics) of formally infinite densities which are ren- 
dered finite mainly by the velocity dispersion of dark mat- 
ter. 

In this paper, although we consider an Einstein-de Sit- 
ter Universe, we do expect our results to give reasonable ap- 
proximation for a ACDM Universe as well. The role of dark 
energy becomes significant at rather small redshifts (~ 0.2) 
which we expect to be well after the formation of the typi- 
cal dark matter haloes we consider here. Furthermore, once 
a particle turns around and collapses, it separates from the 
background expansion and its subsequent motion should not 
be affected by the A term. 

The effect of angular momentum is not considered here 



and we do not expect it to have a significant effect far from 
the core of dark matter haloes. Indeed detailed studies show 
that angular momentu m will not alte r the results of the self- 
similar model (e.g. see (iNussedbOOlh ). 

In this model, the initial perturbation is assumed to be 
scale free, i.e. 



SM 
5 = — — (X M 
M 



at 



t = ti 



(1) 



Such scale-free perturbations result naturally in cosmologi- 
cal models with scale-free initial spectra. For a power-law 
spectrum P{k) ~ k" , the mass variance a{M) scales as 
j(^-(n+3)/6^ which yields 

n-h3 



6 



(2) 



As more and more mass shells turn around and collapse, 
the mass of halo grows at a rate of t^^^'^ . The problem ex- 
hibits similarity solution asymptotically and can be solved 
numerically. The numerical solutions confirm the theoreti- 
cal result s which have been obtained assuming adiabatic in- 
variance (jFillmore fc GoldreichI [T984 ) . Asymptotically, the 
density profile reaches the power-law : 

-2 



P = 



^-9e/(l + 36) 



(3) 



For galactic-scale fluctuations (n ~ 1.2 as given by the 
galaxy-galaxy correlation function) one has e ~ 0.2 — 0.3. 
Here, we set e = 0. 2, also for the purpose of com paring a few 
of our results with iFiflmore fc GoldreichI () 19841 ) . The evolu- 
tion equation, namely Newton's equation of motion written 
in self-similar variables, is0 



1 4 
3^9^ 



dX 
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2 2 
3^9^ 



A 



where 
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A : 



2 /2 
exp — (e 



r 

rta 
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and 



1)C Af(A) 



^ = ln 



tt. 



(4) 



(5) 



The dimensionless radius A is given in terms of the physi- 
cal radius r and present turnaround radius rta and ^ is the 
dimensionless time given in terms of the turnaround time 
tta for each particle. Equation Q has been solved with the 
initial condition given at f = (corresponding to t — tta) 

P , (6) 



A = 1 and — ^ , 
with a prior knowledge of the mass M(A) and where 



P = 



1 + 



3e 



(7) 



The solution to Q can been obtained iteratively : following 
an initial guess for M(A), equation (|4]) is integrated, then 
M(A) is evaluated, and the procedure is continued until the 
desired level of convergence is achieved. Here, we take a sim- 
pler approach. At small values of A (A <C 1), mass becomes a 



1 This is equation (23) of lFillmore fc GoldreichI lll984l ). rewritten 
in variables JS} [also used in iBertschjngeiTllQSsI )]. 
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Figure 1. Particle trajectories in real and phase spaces are shown for two values of the parameter e where A is the non-dimensional radius 
and 5 is the non-dimensional time. The curves for e = 1 are featured in order to highlight the phase transition that occurs at e = 2/3. 
The physical origin for that transition is that for e < 2/3, at late times, the mass within a given radius r is do minated by the contribution 
from particles from outside of r. In contrast for e > 2/3, the mass is dominated by particles within r (see jFillmore &: GoldreichlfigsJ) 
for further details). 




10'^ 10"' 10"^ 10"' 

Figure 2. The left panel features the total numerically-obtained density for e = 0.2 [solid black curve, using expression l(9j]. This 
solution is compared in the right panel to our analytic fit [dashed red curve, using expression l|19| l] . 



power- law M{\) A (Fillmore & Goldreich 1984). We thus 
fit M{\) for e = 0.2 by the following profile 



M{X) 

Mta 



A 



1 + A 



(8) 



where Mta ~ (Sn/A)^. The solutions are shown in Fig. [T]for 
two choices of e. The approximation ([8} generates relatively 
small errors. A notable discrepancy between the results ob- 



tained through approximation ^ and by iteration appears 
only at the relatively large values of A ~ 1. However, for 
all the caustics under consideration the value of A is far less 
than one. The largest value of A for the first caustic is at 
A ~ 0.36 for e = 1 and at A ~ 0.12 for e = 0.2 (see the ta- 
ble in Fig. Al of iMohavaee fc Shandariiil (|2006l ) for further 
details). 
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Figure 3. The dotted lines are the fittings [presented in expressions (|14| l to l|17p ] to numerical results from integration of the differential 
equation The black open circles and the green crosses correspond respectively to e = 0.2 and e = 1. 



The density profile is given by 
8A2 



pb 



(9) 



where pt is the cosmological background density which we 
take to be 1.6 x lO"*^ GeV cm and Us. is the turnaround 
tim e of the particle tha t is at the jth point where A = A(5) 
see iBertschingeil (|l985l ) for further explanation]. 

The density (|9]) is evaluated numerically and plotted 
in Fig. [2] Evidently, the density profile is convolved with 
many spikes, namely caustics, where the density formally 
diverges. However, there is a natural cut-off to the density 
of the caustics due to the finite velocity dispersion of dark 



matter jMohavaee fc Shandaririll2006h . For a generic value 
of e, the thickness of the caustic shell may be expressed as 



(Stt 



,2/3 



1/3 



g(F-l/3)5fc 



(10) 



4 \3e + lj rt: 
The maximum density at the caustic positions is given by 

Gk 



Pk.max 



Pb 



where 



5/6 ^(6-9F)5fc/3 



(11) 



(12) 



32y=2Aj V3e + 1. 
and t is the age of the Universe, a is the present-day velocity 
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k 






K 


Afc 


1 


0.86369 


0.12460 


-8.8 


-0.0222 


2 


1.09899 


0.06987 


-16.2 


-0.0066 


3 


"1 A r\r\r\ 

1.24000 


0.04899 


-23.4 


-0.0032 


4 


1.34200 


0.03782 


-30.5 


-0.0019 


5 


1.42200 


0.03083 


-37.6 


-0.0012 


6 


1.48800 


0.02604 


-44.6 


-0.0009 


7 


1.54499 


0.02254 


-51.6 


-0.0007 


8 


1.59399 


0.01988 


-58.6 


-0.0005 


9 


1.63800 


0.01778 


-65.5 


-0.0004 


10 


1.67700 


0.01608 


-72.5 


-0.0003 



Table 1. The caustic parameters, for e = 0.2, obtained from nu- 
merical solution of l|4]l for the first ten caustics. Analytic fits I I14I I- 
((TtJ reproduce these results and give their values for the rest of 
the inner caustics. 



dispersion of dark matter particles which is that at decou- 
pling re-scaled with the expansion factor and ^fc = \n{t/tta) 
is the dimensionless time given in terms of the turnaround 
time tta of the particle that is now in caustic k. The first ten 
values of these parameters are listed in Table [T] assuming a 
velocity dispersion a of 0.03 cm s"^. 

The density profile of caustics is simply given by 



Pk 



pk. max 



Xk 



lAAfcl ^ A s; Afc 



\A\k\ 
Afc- A 



Pfc,n 



Afe+1 ^ A ^5 Afc — lAAfc 



(13) 

For e — 0.2, the four caustic parameters can be obtained 
for all caustics from the following fits (see also Fig. [3]) ; 



Afc = 0.14315 

^k = 0.8169 + 0.37735 In(fc) , 
a; = = -8.0272 



Afc = 



-0.028533 k' 



-1.9439 



(14) 
(15) 

(16) 
(17) 



Similar fits can be obtained for other values of e. The lower 
envelope (LE) of the density curve (dashed blue curve in 
Fig. [21) is fitted by 

-2.1093 



Ple 

Pb 



= 0.80527 A 



(18) 



Hence the total density curve can be fitted using the above 
expression together with the expressions (|13p for the density 
profile pfc of the peaks, i.e. 

p = Ple + pk . (19) 

In the right panel of Fig. O the above analytic fit is featured 
by the dashed red curve while the true density Q from 
the numerical solution to Q corresponds to the solid black 
curve. The fitted density profile is an excellent description 
of caustics even if it differs slightly from the true density 
profile far from the caustics. 



3 THE COSMIC RAY SIGNAL OF DARK 
MATTER CAUSTICS 

(i) Cosmic ray propagation : the salient features 



Supersymmetry and extra-dimension extensions of the stan- 
dard electroweak model provide a natural candidate for dark 
matter in the form of a weakly interacting and massive parti- 
cle (hereafter WIMP). These species should fill up the galac- 
tic halo. They should also annihilate in pair and produce a 
host of particles among which are rare cosmic rays such as 
antiprotons and positrons. The production rate q (x, E) of 
these antimatter particles depends on their energy E and is 
related to the WIMP annihilation cross section (Tann through 



q (x, E) =r] (a^nnv) 



P(x) 



f{E) 



(20) 



The coefficient 77 is a quantum factor equal to 1/2 for a 
self-conjugate particle like a Majorana fermion or to 1/4 
otherwise. In what follows, we borrow an example from su- 
persymmetry and set 77 = 1/2. The annihilation cross section 
is averaged over the momenta of the incoming DM particles 
to yield (a^^nnv) whose value depends on the specific micro- 
scopic interactions at stake. The WIMP mass is denoted by 
. The energy distribution of the positrons dN^+ / dE^+ 
or of the antiprotons dNp/dEp produced in a single annihi- 
lation is generically denoted by f{E). 

Whatever the source mechanism, charged cosmic rays 
subsequently propagate through the galactic magnetic field 
and are deflected by its irregularities : the Alfven waves. In 
the regime where the magnetic turbulence is strong (which 
is the case for the Milky Way) cosmic ray transport needs 
to be investigated numerically. Monte Carlo simulations 
jCasse et al.ll200"3 ) indicate that it is similar to space dif- 
fusion with a coefficient 



K{E) =Ko 13 71^ 



(21) 



which increases as a power law with the momentum to elec- 
tric charge ratio TZ = p/q - also called the rigidity - of the 
particle. In addition, because the scattering centres drift in- 
side the Milky Way with a velocity Va ~ 20 to 100 km s~^ , a 
second order Fermi mechanism is responsible for some mild 
diffusive re-acceleration. Its coefficient Kee depends on the 
particle velocity /3 and total energy E and is related to the 
space diffusion coefficient K{E) through 



Kee ~ 



2 2 

9 " K{E) 



(22) 



In the case of positrons, diffusive re-acceleration is com- 
pletely dominated by energy losses. Finally, galactic convec- 
tion wipes cosmic rays away from the disc with a velocity 
Vc ~ 5 to 15 km s~\ 

After this short digest of cosmic ray transport, we can 
write the master equation fulfilled by the space and energy 
distribution function -0 = dn/dE as 

{Vc Tp) - KAiP + (23) 
+ Oe [b'°''{E) ^ - Kee{E) dEi>] = q (x, E) . 

This equation applies to any species (protons, antiprotons 
or positrons) as long as the rates for production q and en- 
ergy loss b^°'"'{E) are properly accounted for. It has been 
solved within the framework of the semi-analytical two- 
zone model which has been ext e nsively discussed in pre- 
vious works (|Maurin et al.l I2OOII : iDonato et"al] bOOlh and 
whose salient features we briefiy recall now. According to 
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Case 


5 


Ko [kpcVMyr] 


L [kpc] 


Vc [km/s] 


Va [km/s] 


MIN 


0.85 


0.0016 


1 


13.5 


22.4 


MED 


0.70 


0.0112 


4 


12 


52.9 


MAX 


0.46 


0.0765 


15 


5 


117.6 



Table 2. Typical c ombinations of diffusion parameters that are compatible with the B/C analysis llMaurin et al.ll200lh . As shown in 
l|Donato et al.il2004h . these propagation models correspond respectively to minimal, medium and maximal primary antiproton fluxes. 



our approach, a steady state is assumed and the region of the 
Galaxy inside which cosmic rays diffuse (the so-called diffu- 
sive halo or DH) is pictured as a thick disc which matches 
the circular structure of the Milk Way [as shown in Fig. 2 
of jMaurin et al.ll2002l 'l]. The galactic disc of stars and gas, 
where primary cosmic rays are accelerated, lies in the mid- 
dle. It extends radially 20 kpc from the centre and has a 
half-thickness h of 100 pc. Confinement layers where cosmic 
rays are trapped by diffusion lie above and beneath this thin 
disc of gas. The intergalactic medium starts at the vertical 
boundaries z = ±L as well as beyond a radius of r = i? = 20 
kpc. Notice that the half-thickness L of the diffusive halo is 
not known and reasonable values range from 1 to 15 kpc. 
The diffusion coefficient K is the same everywhere whereas 
the convective velocity is exclusively vertical with compo- 
nent Vciz) = Vc sign(z). The galactic wind, which is pro- 
duced by the bulk of the disc stars like the Sun, drifts away 
from its progenitors along the vertical directions, hence the 
particular form assumed here for Vc- Notice also that the 
normalization coefficient Ko, the index S, the galactic drift 
velocity Vc and the Alfven velocity Va are all unknown. 
This situation may be remedied with the help of the boron 
to carbon B/C ratio which is quite sensitive to cosmic ray 
transport and which may be used as a constraint. The three 
prop agation model s featu red in Table [2] have been borrowed 
from lDonato et all (|2004l '). The MED configuration provides 
the best fit to the B/C measurements whereas the MIN and 
MAX models lead respectively to the minimal and maximal 
allowed antiproton fluxes which can be produced by WIMP 
annihilation. 



factored out in 



The solution of the master equation p3p may be gener- 
ically expressed as the integral 



{Q,E) = [dEsf (f yis G i^Q, E ^ xs, Es) q{^s,Es). 

J J DH 



(24) 

The energy Es at the source runs over a range which de- 
pends on the nature of the cosmic ray species as discussed 
below. The space integral is performed over the diffusive 
halo. The convolution (|24|l involves the Green function G 
which describes the probability for a cosmic ray that is pro- 
duced at location xs with the energy Es to reach the Earth 
where it is detected with the degraded energy E. The cosmic 
ray space and energy density tp can be translated into the 
differential flux $ = Ptp/i-K where P stands for the particle 
velocity. This flux is expressed in units of particles m~^ s~^ 
sr~^ GeV~^. In the case of WIMP annihilations, it can be 
recast as 



^{Q,E)^T jdEsf{Es)I{E,Es) , 



(25) 



(26) 



The energy distribution f{Es) describes the source spec- 
trum and depends on the WIMP properties. The informa- 
tion on the galactic DM density profile p as well as on the 
propagation of cosmic rays within the Milky Way is summa- 
rized in the halo integral 



I{E,Es)= /d='xsG(xQ,S^xs,Ss) 

J DH 



P0 



(27) 

where the solar neighbourhood DM density is denoted by 
Pq. The halo integral I{E,Es) plays a crucial role in deriv- 
ing the fiux at the Earth of the antimatter species produced 
inside the galactic DM halo by WIMP annihilations. The 
reach of the Green function G depends on the nature of 
the cosmic ray particle (either positrons or antiprotons in 
our case) as well as on the energies E and Es- It delineates 
the region of the Milky Way from which most of the signal 
detected at the Earth originates. This so-called horizon, be- 
yond which the Green function is suppressed, plays a crucial 
role in our discussion as it may or may not reach down to the 
centre of the Milky Way and probe its dense DM content. 

Except for high-energy positrons, the horizon is much 
larger than the separation between adjacent caustics for 
which we find a typical distance of ~ 130 pc in the solar 
neighbourhood. If the DM halo of our Galaxy results from 
the self-similar accretion which we have discussed in Sec- 
tion[2l the WIMP density that enters in the DH integral ((27)l 
is given by the sum of the lower envelope density (|18|l and 
the caustic density (|13p . The latter varies by several orders 
of magnitude at the crossing of each shell as shown in Fig. (2] 
Because the typical scale over which the Green function G 
changes appreciably is much larger than the shell thickness 
Arfc and the shell separation — ''fe+i, we may actually 
average the density pk of the k**^ caustic over the distance 
extending from the inner radius r^+i to the outer boundary 
rfc. This leads respectively to the linear average 



pk ~ pk,msx {2y/ek — eu) , 

as well as to the quadratic average 

{efc(l-lnefe)} 



1? - r? 
Pk — Pfc.max 



(28) 



(29) 



where ~ AAfe/ (Afc — Xk+i)- The DM density squared that 
enters in the halo integral I{E, Es) may be expressed as the 
average 



where the information related to particle physics has been 



p2 = p^E + 2pLEP^ + pI - (30) 

The corresponding root mean square density is featured in 
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1.95 - ^ 




root mean square 
total 

lower — envelop 



r [kpc] 

Figure 4. Typical DM density in the self-similar accretion model 
of Section |2] The solar system is located at = 8 kpc from the 
galactic centre and is surrounded by two nearby thin shells. The 
long-dashed dotted red curve features the root mean square den- 
sity used in the calculation of the halo integral I(E, Eg) whereas 
the short-dashed blue line corresponds to the lower envelope p\_^. 



Fig. |4] as a function of the galactocentric distance r (long- 
dashed dotted red curve) together with the actual DM den- 
sity Ple + Pk (solid black line) and the short-dashed blue 
lower envelope ple- The turnaround radius rta has been de- 
termined by requiring that the linearly averaged DM den- 
sity p = Ple + pfc is equal to a solar neighbourhood value of 
P0 = 0.3 GeV cm~'^. We infer a turnaround radius rta of 
2.7658 Mpc. The solar system is located between the SS**^ 
and 59"^ shells whose inner densities reach a maximum value 
p max — 2.13 X 10^ pb of 34.2 GeV cm . These caustics are 
extremely thin though, and we find AA58/(A58 — A59) ~ 
1.5 X 10"* with a shell separation of 131 pc and a shell 
thickness of AA58 = 7 x 10"" 



The DM shell density derived from relation (|30|l is very 
close to the lower envelope density ple. Both diverge at the 
galactic centre with a profile index of 7 ~ 2.11. This is 
physically unacceptable. The rotation velocity decreases in 
the bulge as the galactocentric distance goes to zero. The 
self-similar model is not a good description of the Milky 
Way DM halo close to its centre and a cut-off radius rcut-off 
needs to be imposed by hand in order to get a more accept- 
able behaviour in that region. The existence and the value 
and origin of a core at the centre of dark matter haloes 
remains an open problem. The core radius claimed in the 
literature spans many orders of magnitude and hence re- 
main uncertain. It depends on very many different parame- 
ters : baryons, angular momentum, presence or absence of a 
central black hole, DM self-annihilation, DM velocity disper- 
sion, merger and tidal efi^ects. We thus present our results for 
a wide range of core values. The maximum range of 3 kpc 
chosen here is motivated by recent studies such as those 
on the core radii of dwarf galaxies (see e.g. (jCoerdt et al.l 
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Navarro, Frenk & White 
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(Navarro et al. 1997) 



Table 3. Parameters for different dark matter density profiles 
[equation ([STJ] for the Milky Way. 



l2006h l. Because baryons dominate the inner potential well 
of the Milky Way, the self-similarity of the DM distribution 
breaks down inside the bulge. A plausible value of ~ 1 kpc 
for the core radius r-cut-off is then set by the bulge extension. 
The case of a 300 pc cut-off radius has also been investigated 
though it may represent a somewhat extreme situation. 

(ii) The positron signal 

The WIMPs which annihilate inside the thin and highly con- 
centrated DM shells of the Milky Way halo are expected to 
yield an a priori much higher positron or antiproton fiux at 
the Earth than in the more conventional case of a smooth 
DM galactic distribution. Previous investigations have con- 
centrated essentially on the models listed in Table [3] for 
which the DM profile is generically given by 



1 



(r/r 



The galactocentric distance r© of the solar system has 
been set equal to 8 kpc. Th e isothermal cored distribu- 
tion (|Bahcall fc Soneira|[l980l ) features a constant DM den- 
sity within th e inner 4 kpc wherea s the Navarro, Frenk and 
White profile [Navarro et al.l (|l997l ) hereafter NFW] diverges 
like at the galactic cen tre. Other numerical investiga- 
tions ([Diemand et al.l |2004| ) find a DM profile index 7 of 
1.16 ± 0.14 in good agreement with the NFW result and 
should lead to the same positron or antiproton flux at the 
Earth. In order to compare the positron signal ^shcu (0, E) 
from DM caustics to the flux "I'smooth (0, E) yielded by the 
DM distributions of Table O we have computed their ra- 
tio. The positron and antiproton signals from annihilating 
WIMPs have already been thoroughly investigated in the 
literature and shown to depend on both the assumed DM 
distribution and the speciflc features of the particle physics 
model selected to describe the DM species. Our aim is here 
to gauge solely the influence of the DM shells on the anti- 
matter cosmic ray signal and to investigate whether or not 
the annihilation of WIMPs is enhanced inside the caustics 
and leads to a larger signal than in the case of an isothermal 
cored proflle or a NFW cusp. 

Because the flux depends on the source spectrum 
f{Es), we have decided to concentrate on the halo integral 
I{E, Es) which specifically encodes information on both the 
DM distribution and cosmic ray propagation. Positrons es- 
sentially diffuse on the irregularities of the magnetic field 
and lose energy through synchrotron radiation and inverse 
Compton scattering on the cosmic microwave background 
radiation and on the galactic starlight. The energy loss rate 
^lossj-^^ increases with the energy E as 

2 



6'°==(£;) 



where 



Eoe'/TE , (32) 
E/Eo and Eq = 1 GeV. The typical energy loss 
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Figure 5. Positron boost factor in the presence of DM caustics normalized to a NFW profile (solid) or to an isothermal cored distribution 
(long-dashed) in the case of the propagation model MED of Table [2] In the left panel, a positron line at Eg = 1 TeV has been assumed 
and is plotted with respect to the energy E at the Earth. In the right panel, the boost factor is featured as a function of the diffusion 
length Ad- Three different values for the cut-off radius of the DM shell distribution have been considered. 



timescale is rg = 10^^ s and m ay be combined (jLavalle et al.l 
l2007al : [PeTahave et al.l \200t ) with the diffusion coefficient 
K{E) to yield the typical diffusion length 

(33) 



The parameter Ko is the normalization of the diffusion coef- 
ficient while f = t{E) — t{Es) is the typical time (including 
both energy losses and diffusion) during whi ch the positron 
energy decreases from Es to E. As shown bv lBaltz fc Edsia 
(|l999h . the positron energy E may be translated into the 
pseudo-time 



i{E) 



TE 



(34) 



where 5 denotes the spectral index of the diffusion coefficient 
K - see relation (|2ip . The diffusion length Ad measures ac- 
tually the extension of the positron sphere. It gauges how far 
positrons travel before being detected at the Earth. A rapid 
inspection of equation (|34|l indicates that Ad increases as 
the detected energy E decreases except for energies Es at 
the source very close to E. It is well known that the positron 
sphere is fairly reduced at high energies, say above ~ 100 
GeV, whereas it extends over several kiloparsecs below 10 
GeV. 

The transition occurs between 10 and 100 GeV as fea- 
tured in the left panel of Fig. [5] where the positron boost, 
which we define as 



4'hcll(-£', Es)/Isino 



th{E, Es 



(35) 



is plotted as a function of the positron energy E at the 
Earth {i.e. after propagation and energy loss) for a fixed 
value of the injection energy Es- Such a situation arises for 
Kaluza-Klein species because they can annihilate directly 
into electron-positron pairs. The injection energy Es is equal 
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Figure 6. Same as in the left panel of Fig.[5]but now for a cut-off 
radius of 300 pc and for the three propagation models of Table [2] 



to the WIMP mass and a value of 1 TeV has been chosen 
here, in agreement with universal extra dimension (UED) 
models. The solid and long-dashed curves correspond re- 
spectively to the NFW profile and isothermal cored distri- 
bution of Table |3] Three values of the core radius rcut-off 
are displayed. For energies at the Earth exceeding ~ 100 
GeV, the positron sphere does not extend far away. The sig- 
nal originates from the solar neighbourhood where the DM 
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density is equal to p©. The halo integrals P {E,Es) are 
identical whatever the underlying model for the DM density 
profile. As expected, the positron boost factor is approxi- 
mately equal to unity. Below ~ 100 GeV, the positron sphere 
widens as E decreases and Ad consequently gets larger. Be- 
cause DM is more abundant in the shell scenario than what 
is expected in the case of the smooth profiles of Table |3] the 
positron boost increases as the energy E goes down from 
100 GeV to 1 GeV. This trend is particularly pronounced 
for the isothermal sphere (long-dashed curves) for which the 
DM abundance is the smallest in the inner parts of the Milky 
Way. The boost factor is shifted downwards (solid lines) by 
a factor of ~ 2 when a NFW profile is considered instead. As 
the cut-off radius is decreased from 3 kpc down to 300 pc, 
the shell DM density gets larger at the galactic centre and 
so does the boost. Notice also that the true argument of the 
halo integrals P is the diffusion length Ad which combines 
both energies E and Es- The right panel of Fig. [5l where 
the positron boost is presented as a function of Ad, is a 
mere reflection of the left plot. Because the diffusion length 
increases as the energy E at the Earth goes down, these 
panels are indeed symmetrical. 

The solid and long-dashed red curves of Fig. [5] corre- 
spond to a cut-off radius of 300 pc. They have been plotted 
in Fig. |S]as a function of the positron energy E at the Earth 
for the three different propagation models of Tabled In the 
MAX configuration, the spectral index 5 is the smallest and 
the pseudo-time t{E) is most sensitive to the energy E as 
expected from relation (|34p . This trend is strengthened by 
a large diffusion coefficient Kq. The positron horizon probes 
the galactic centre as soon as the energy E at the Earth gets 
lower than ~ 50 GeV. The same regime is reached for an en- 
ergy of ~ 10 GeV in the case of the MED model whereas the 
positron sphere never extends to the galactic centre in the 
MIN case even for energies at the Earth as low as 1 GeV. 

The positron line illustrates the behaviour of the hori- 
zon as a function of the energy. This case has been dis- 
cussed for pedagogical purposes though it is inspired by 
Kaluza- Klein models where WIMPs annihilate directly into 
electron-positron pairs. The boost can actually reach large 
values as -Es is kept fixed whereas the detection energy E 
is decreased down to a few GeV. We should also expect a 
continuous positron spectrum to arise from WIMP annihi- 
lations. Its distribution f{Es) decreases more or less mildly 
with the injection energy Es - The signal at the Earth is then 
given by the convolution of f{Es) with the positron Green 
function. The positron sphere is on average smaller than for 
the line and we definitely expect smaller boost factors whose 
actual values depend on the steepness of the injection spec- 
trum. 

(iii) The antiproton signature 

The propagation of cosmic ray antiproton is dominated by 
diffusion. Energy losses as well as diffusive re-acceleration do 
not play any major role. A very crude approximation for the 
antiproton Green function is obtained by neglecting galac- 
tic convection and solving the resulting Poisson equation in 
infinite space. This yields the antiproton propagator 
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(36) 



A-kK{E) r® ' 

where rm denotes the distance between the Earth and the 
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Figure 7. Antiproton boost factor Bp in the presence of DM 
caustics with respect to a NFW profile. A cut-ofF radius of 300 pc 
has been assumed for the shell distribution. The three propaga- 
tion models of Table [2] are featured. The boost factor is plotted 
as a function of the interstellar antiproton energy Tp^ . 



source. The sole merit of this expression is to exhibit the 
importance of remote sources. We will therefore keep in 
mind that the antiproton sphere is more extended than 
for positrons. The finite thickness of the diffusive halo is 
nevertheless a limiting factor since cosmic rays may escape 
through the vertical boundaries as they wander towards the 
Earth. As a consequence, the size of the antiproton sphere 
cannot be much larger than the DH half-thickness L. Galac- 
tic convection comes also into play. If the wind velocity 
Vc is large, cosmic rays are efficiently blown outside the 
Milky Way. This process limits further the reach of the an- 
tiproton sphere. The solution of the master equation (|23|l 
has bee n thoroughly invest i gated and several different tech- 
niques (jMaurin et al.ll2006l : iBringmann fc Salatill2007t l lead 
essentially to the same fluxes at the Earth. In the Green 
function formalism, the diffusive halo is pictured as an in- 
finite slab with no radial boundaries. On the contrary, the 
Bessel expansion method takes advantage of the axial sym- 
metry of the propagation region and enforces a vanishing 
cosmic ray fiux at a distance _R = 20 kpc from the rotation 
axis of the Milky Way. 

Since energy losses and diffusive re-acceleration are neg- 
ligible, the antiproton energy Es at the source is equal to its 
energy E at the Earth. Expression (|25|l may be considerably 
simplified into 



<I>p {Q,E)=Tf{E)IiE) 



(37) 



where the halo integral I is now a function of the antiproton 
energy E = Es- The signal generated by the DM shells of 
Section [2] is enhanced with respect to the smooth distribu- 
tions of Table |3] by the boost factor 



4on(S)/Cooth(S) 



(38) 
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Figure 8. Antiproton boost factor in the presence of DM caustics normalized to a NFW profile (solid) or an isothermal cored distribution 
(long-dashed) in the case of the propagation model MED of Table [2] Three diff'cront values for the cut-ofF radius of the DM shell 
distribution have been considered. Energy losses and diffusive re-acceleration have been taken into account in the right panel as well as 
the production of tertiary antiprotons which induces a dramatic increase of the boost at low energy. 




Figure 10. Signals expected in the case of a 1.7 TeV Wino and a DM shell distribution with a cut-off radius of 300 pc. In the left panel, 
the antiprot on flu x t&p is plott ed as a function of the antiproton kinetic energy tTOA ^^^^ compared to several measurements 
jOrito et al.ll2000l ; iMaeno et ai]|200ll : [Boezio et al.ll200ll : IXguilar et al.ll2002l ) and projected observations jPicozza et al.ll200l iBaraol 
l2004h . The MED propagation model has been assumed. In the righ t panel, the positron fraction is featured for the MAX and MED 
propagation sets together with the HEAT data llBarwick et al.|[l997l ). 



The latter is plotted in Fig. [7] as a function of the antipro- 
ton interstellar kinetic energy for the three propagation 
models of Table [2] A cut-ofF radius of 300 pc has been as- 
sumed for the DM shell distribution. Its boost Bp is cal- 
culated with respect to a smooth NFW halo. The Green 



function technique and the Bessel expansion method yield 
identical results. Energy losses and diffusive re-acceleration 
have been switched off as well as the production of tertiary 
antiprotons. The MAX propagation model (red curve) is 
characterized by a large value of the half-thickness L and a 
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Figure 9. Same plot as the left panel of Fig. [8] where the prop- 
agation model is now the MAX configuration of Table [2] 



small convection velocity Vc ■ The antiproton sphere spreads 
therefore over a large portion of the Milky Way and reaches 
the inner dense regions of its DM halo. The boost factor 
is essentially constant and is equal to ~ 30 above a few 
GeV. For exactly the opposite reasons, the MIN configura- 
tion (blue line) does not lead to any enhancement of the 
antiproton signal in the presence of DM caustics. The dif- 
fusive halo is too thin and the galactic wind too strong to 
let the antiprotons originating from the galactic centre to 
reach the Earth. The MED case (in green) features the in- 
termediate situation. Above ~ 100 GeV, diffusion takes over 
galactic convection and the antiproton sphere is essentially 
limited by the half-thickness L of the DH. That sphere ex- 
tends sufficiently close to the centre to become sensitive to 
its large DM density, hence a boost factor identical to the 
MAX result. Below ~ 10 GeV, the galactic wind is on the 
contrary strong enough to restrain the antiproton sphere 
from spreading and the MIN situation is recovered. 

In Fig.m the MED propagation mode has been selected. 
The boost Bp is featured as a function of the antiproton in- 
terstellar kinetic energy for three different values of the 
cut-off radius rcut-ofi- It has been normalized to a NFW pro- 
file (solid) or an isothermal cored distribution (long-dashed). 
The solid red curve of the left panel corresponds therefore 
to the green line of Fig. [7] and the decrease of the boost at 
low energy can be explained as a result of the galactic wind 
taking over diffusion. In the right panel, the Bessel expan- 
sion technique alone has been used in order to implement 
energy losses and diffusive re-acceleration. The inelastic and 
yet non- annihilating interaction of an antiproton with an hy- 
drogen atom of the interstellar medium is possible through 
the production of a A resonance which subsequently decays 
into a proton and a pion. The antiproton is not destroyed 
but loses a significant portion of its energy in a single colli- 
sion. This so-called tertiary process has also been taken into 



account. The curves are significantly flattened since the low- 
energy tail of the spectrum is replenished by higher energy 
antiprotons undergoing these inelastic and non- annihilating 
scatterings. Since the energy E drifts now with respect to 
its initial value Es, the injected spectrum matters. In the 
right panel, f{Es) has been assumed to be constant. This 
is an unrealistic situation as the antiproton spectrum gen- 
erated by WIMP annihilation is expected to be softer than 
a flat distribution. The panels of Fig. [8] feature therefore 
two extreme trends and the actual boost should follow an 
intermediate behaviour. 

In Fig.[9l energy losses and diffusive re-acceleration have 
been switched off as well as the above mentioned tertiary 
process. The situation is therefore similar to what has been 
assumed in the left panel of Fig. [8] up to the propagation 
model which is now the MAX configuration of Table (2] The 
right panel of Fig. |S] (MED propagation with tertiaries and 
diffusive re-acceleration) is similar to Fig. [9] (MAX propa- 
gation without tertiaries nor diffusive re-acceleration). We 
conclude that in the MED configuration, the antiproton 
horizon may be somewhat stretched at low energy by the 
tertiary process which allows a fraction of the high energy 
antiprotons to make it from the inner Galaxy to the Earth 
where they subsequently populate the energy band below ~ 
10 GeV. 

The annihilation cross section of WIMPs is in gen- 
eral set by the requirement that their relic abundance to- 
day should be compatible with the WMAP measurement of 
Q.M ~ 0.25. This leads to a value of {(Tanni') of order a few 
lO"'^^ cm^ s~^ and to an undetectable signal. The presence 
of substructures in the Milky Way halo has been invoked to 
increase the cosmic r ay fluxes at the Eart h by an amount 
which is still debated (|Lavalle et al.|[2007lj ). The aim of this 
article is actually to estimate that value in the case of caus- 
tics. We find that the boost hardly exceeds ~ 30. With such 
a value and assuming the above mentioned annihilation cross 
section, the antiproton and positron signals are barely de- 
tectable. There is however a case where {a^^nnv) could be as 
large as 10~^* cm^ s~^. In an e xtensive analysis of a fair! 
generic supersymmetric model, IProfumo fc Yagunal (|200' 
have shown that a heavy Wino-like LSP often appears as a 
plausible DM candidate. This species is expected for exam- 
ple in ano maly media ted supersymmetry breaking (AMSB) 
scenarios (|Ullid [200 ll ) and relic density requirements point 
towards a mass of 1.7 TeV. We emphasize that this config- 
uration is quite generic. Non-perturbative, binding energy 
effects imply a greatly enhanced annihilation cross section 
today, when the Winos ha v e very small galactic velocities 
as shown by iHisano et al.l (|2005r ). In Fig. \M a value of 
1.02 X 10~^^ cm^ s~^ has been assumed for {(Jann'j)- In this 
limit, heavy Winos annihilate almost exclusively into gauge 
bosons, with a branching ratio of 80% into W^W' pairs 
and of 20% into Z^Z^^ pairs in our example. In the left 
panel, the antiproton spectrum is plotted as a function of 
the antiproton kinetic energy. The MED propagation mode 
is selected here and since in that case the antiproton hori- 
zon is fairly limited at low energy, the yield from DM an- 
nihilation is small and compatible with the measurements. 
Above ~ 10 GeV, the signal however increases dramatically 
with respect to the background (yellow band) which con- 
sists of secondary antiprotons. The latter are produced con- 
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ventionally by the spallations of high-energy CR nuclei on 
the interstellar gas. The s olid red curve is ev en compatible 
with the CAPRICE data (|Boezio et alj 120011 ) which possi- 
bly hint at a spectral distortion. The largest discrepancy 
between the DM shell signal and the secondary background 
is reached for an energy of 300 GeV. The antiproton sphere 
clearly reaches the core of the DM shell distribution. Notice 
that no artificial enhancement of the antiproton primary 
flux is necessary here. We have used solely the DM distri- 
bution resulting from the self-similar infall model of Sec- 
tion [5] with no extra boost factor. Our only assumption is a 
cut-off radius of 0.3 kpc. Should we have considered larger 
core radii, we would have obtained however a smaller signal. 
In the right panel, the positron fraction is presented as a 
function of the positron energy E at the Earth. The solid 
magenta curve corresponds to t he ba ckground alone as de- 
rived by iMoskalenko fc Stronj (|l998l ). Secondary positrons 
are produced in the Galaxy from the collisions of cosmic 
ray proton and helium nuclei on the interstellar medium. 
The short-dashed green line features the signal of our 1.7 
TeV Wino and is based on the MED propagation model. 
The positron sphere is sensitive to the high DM density of 
the central caustics only at low energy, in a regime where 
the secondary background is large enough to hide easily a 
spectral distortion. This is no longer the case if the MAX 
propagation model is substituted for the MED configura- 
tion. The positron sphere becomes large enough to reach 
the galactic centre below 50 GeV and the long-dashed red 
curve exhibits a clear exces s with respect to the background. 
The HEAT measurements (|Barwick et al.|[l997l 'l are not well 
reproduced though. They point towards a positron excess 
which increases above 10 GeV, a trend which our shell den- 
sity (|30|l clearly does not match. A Kaluza-Klein DM species 
with a substantial positron line has been invoked to fit th e 
HEAT distortion (|Cheng et al.ll2002l : [Hooper fc Silkll2005h . 
However, an artificial boost factor of 5 x (0.43/0.3)^ ~ 10 on 
the entire positron energy range is necessary. Another possi- 
bility relies on the presence of high concentrations of WIMPs 
in the solar neighbourhood. The increase of the HEAT ex- 
cess with energy indicates that a local mechanism is at play. 
The probability fo r a DM clump to be close to the Earth ha s 
been investigated (|Hooper et al.ll2004l : iLavalle et allbOOTal ). 
In our case, the caustics themselves could play the same role 
as the clumps, because the DM density reaches there a priori 
extreme values. We may wonder then how large the positron 
signal could be, should the solar system be close to one of 
these shells. 



4 DARK MATTER DRAPERIES AND THE 
POSITRON EXCESS 

Unlike antiprotons which have a very long diffusion length, 
the positrons that are detected at the Earth before signif- 
icantly losing their energy [E ~ Es) have covered on av- 
erage a small distance. For this reason, they can be impor- 
tant tracers of nearby density peak structures such as the 
dark matter "draperies" on which we focus here. As men- 
tioned in Section [31 the solar system is located between two 
nearby shells. In our self-similar infall model, the 58"^ and 
SO*'' caustics lie respectively at a distance of 8.0348 and 
7.9040 kpc from the galactic centre. The shell separation 



e = r58 — rsg is equal to 130.8 pc. It is 8 orders of magni- 
tude larger than the shell thickness for which a typical value 
is given by Arsg — 5.9 x 10^ km, hence the term "drap- 
ery". A turnaround radius rta of 2.7658 Mpc has been ob- 
tained by setting the linearly-averaged density p — piE + psi 
equal to the solar neighbourhood value pQ — 0.3 GeV cm~"^. 
The density psg of the 58"^ caustic reaches a maximum of 
P58,max ~ 34.2 GeV cm^'' inside an homogeneous slab whose 
thickness Arsg is equal to the distance between the Sun and 
Mercury. Should the Earth wander in this region of highly 
concentrated DM, the positron signal would be enhanced by 
a factor of (pss.max/p©)'^ ~ 13, 000, hence the importance of 
exploring this possibility in greater detail. 

To start, we note that changing slightly the value of 
P0 modifies the turnaround radius and affects the relative 
position of the solar system and the DM shells. The Earth 
may lie within the densest region of a caustic but the prob- 
ability for such a situation is grossly given by the ratio 
£58 — Ar58/e ~ 1.5 X 10^** and is vanishingly small. 

Playing the devil's advocate, let us nevertheless assume 
that the solar system is indeed embedded inside this density 
peak region with typical DM density psg.max- Should this be 
the case, the positron signal would still be so smeared by 
diffusion and energy losses that it would be marginally en- 
hanced with respect to the results presented in Section |3l 
Even in the extreme situation of a line where the injec- 
tion energy is Es = m^, positrons are detected at energies 
E < Es and originate from a sphere all the more extended 
as the energy difference AE = Es — E is large. Let us con- 
sider the MED propagation model. The spectral index of the 
diffusion coefficient K is 5 — 0.7. In the limit where the en- 
ergy difference AE is small compared to the positron energy 
E ~ Es, we may differentiate relation (|34|l and combine the 
result with definition (|33|l to get 

Ad = 3.77 kpc Ae'/^ eg° '^'^ , (39) 

where es ~ Es/ Eq and Eo = l GeV. For a typical positron 
energy of 100 GeV, the diffusion length is equal to 



Ad = 1.89 kpc ^— , (40) 

Because the experimental determination of the energy E 
is performed with a limited accuracy, positrons emitted 
monochromatically at Es = 100 GeV are collected in an en- 
ergy band set by the resolution of the instrument. Even in 
the optimistic case where energy is measured at the 1% level. 
I.e. with AE/Es = 10^^, relation ([301 yields a positron 
sphere which reaches as far as 190 pc from the Earth at E 
= 99 GeV. The signal collected in the energy bin of the line 
does not originate solely from the densest part of the caus- 
tic. It has been produced from a much more extended zone 
whose average DM density is 

p ~ PLE ^ P0 = 0.3 GeV cm"^ . (41) 

Remember that if the shells are dense, they are also ex- 
tremely thin. The DM density p58 for instance undergoes a 
discontinuity at the external boundary of the 58**^ caustic. It 
also drops violently inwards, decreasing from 34.2 down to 
0.3 GeV cm^'' on a distance of 2.5 x 10^^ pc only. Averaged 
over a distance D, the 58**^ caustic density is given by 

P58 = 2 P58,max \ , (42) 
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where relation ((28} has been used. For a distance D — 1 pc, 
we get a value of 0.095 GeV cm"'^ which barely represents a 
third of the solar neighbourhood value. Averaged now over 
the distance e separating two nearby shells, it amounts to 
2.8% of the density p©, in agreement with the relative po- 
sitions of the long-dashed dotted red curve (averaged total 
DM density) and the short-dashed blue line (lower envelope) 
of Fig. H 

In order to elaborate on the smearing of the positron 
horizon by the limited resolution of energy measurements, 
some modelling is necessary. Since the shells are extremely 
thin, we will replace the actual DM density by the fine- 
grained distribution peff where 



and 



2 2 

-2- = -2- + 2^ AeS{z- 

"O "O shrlls fc 



Zk) 



(43) 



The shells are now pictured as infinitely thin slabs located at 
the positions z = Zk on the radial axis connecting the galac- 
tic centre to the Earth. They are separated from each other 
by the distance e defined above. Setting the total number of 
DM annihilations constant, the parameter A is found to be 
equal to 



A : 



Pi 



P5S 
PQ 



(44) 



with a typical value of 0.0554. The effect which we dis- 
cuss here is the strongest in the case of a positron line for 
which the injection spectrum is f{Es) ~ S {Es — m^). The 
positron flux at the Earth is given by the product 



<E>,+ (0,i?)=^r {E,Es=m^) 



(45) 



where the halo integral J*^ is the convolution (|27|) between 
the DM density squared and the positron Green function 



. Close to the line energy, for E 



the positron 



horizon is so limited that we may not replace the actual 
DM density by its coarse-grained average (|30|l . especially if 
the Earth lies in the region where the caustic density is the 
highest. The results presented in Section|3]should potentially 
vary if we compute now the halo integral I'' with the fine- 
grained density (|43p . The relative change in the positron fiux 
is equal to the ratio 



sr 



^ ^coarse 



Because the positron propagator may be expressed as 



(x,£; ^ xs,£'s) 



Eoe- 



G (x <— xs; Ad) 



(46) 



(47) 



the halo integral J*^ is related to the convolution / of the 
DM density squared with the heat Green function G through 



r {E,Es) 



TE 

Eoe^ 



/(Ad) 



where 



/(Ad) 



d' xs G (x0 <— xs; Ad) 

DH 



P(xs) 

P© 



(48) 



(49) 



Considering the fine-grained DM distribution (|43p instead 
of the coarse-grained average (|30[l induces a relative change 
in the positron fiux equal to (5///coarse where 



51 = /fine — Ic. 



(50) 



coarse 



= /d='xs(5(xo «-xs;Ad) 4- ■ (51) 

J DH P© 

We are interested here in the limit where the positron diffu- 
sion length Ad is small compared to the intershell separation 
e ~ 131 pc. In this regime. Ad is much smaller than the 
DH half- thickness L and the heat Green function G simpli- 
fies into (|Baltz fc Edsio|[l999l : iLavaUe et al.ll2007ah 



G (x <— xs; Ad) 



.w,. 3 exp. 
tt'^/^Ad L ^d" 



(x - xs) 



(52) 
~ 1. 



Since ~ Pq , we are led to the conclusion that /coai 
Our task consists then in evaluating the difference 

5/(Ad) = /d'xs G(x0 ^xs;Ad) l^^^-^j , (53) 

^ DH I P© J 

and in comparing it to unity. With our fine-grained expres- 
sion (|43p . this integral becomes 

5/(Ad)= fd^xsGA\-l+ eS{z-Zk)\ . (54) 

To proceed further, we simplify the Gaussian expression (|52[) 
for the heat propagator G into the step function 



G = — e (Ad - rg 
Vd 



(55) 



where denotes the distance of the source to the Earth. 
The radius of the positron sphere is Ad and its volume is 
Vd = 47rAD^/3. We readily find that 



<5/(Ad) 



E 



Vd 



(Ad' 



Zk) 



(56) 



where the sum runs over the DM shells that intersect the 
positron sphere. In the limit where the positron diffusion 
length Ad is much larger than the caustic separation e, 
the difference 51 vanishes as expected. In this regime, the 
positron sphere is so large that the coarse-grained den- 
sity (|30p is enough to describe appropriately the DM distri- 
bution. We are actually interested here in the opposite limit 
where the ratio Ad /e is small and where the draperies come 
into play. The difference 51 may be conveniently expressed 



5I{Xu) =A — V, 
Ad 

where the drapery function D is defined as 

^ Ad V- 3 /, 2^ 

e , f^.. 4 V Ad' 



(57) 



(58) 



The latter is plotted in Fig. [TT] as a function of the ra- 
tio Ad/c. The minimal value of the difference SI is reached 
when no DM shell intersects the positron sphere. The Earth 
lies between two nearby draperies. This case is featured for 
instance by the short-dashed green and long-dashed dotted 
blue curves of Fig. [TT] and provided that the positron dif- 
fusion length Ad is small enough. The minimal value of 51 
is equal to — ^ and is negligible with respect to /coarse — 1- 
More complicated is the maximal configuration for which 
the Earth lies exactly inside a drapery and where the DM 
density is the highest. This case corresponds to the solid 
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Figure 11. The drapery function D is plotted against the ratio 
Ad/g for three different values of the distance to the nearest DM 
shell. The solid red curve features the case where the Earth is 
located exactly inside a drapery and where the DM density is the 
highest. The integral &I is maximal. The long-dashed dotted blue 
line corresponds to the opposite situation where the Earth lies 
at mid-distance between two caustics. The integral 51 is minimal 
with a negative value. The intermediate situation for which the 
nearest shell is located at a distance of e/4 leads to the short- 
dashed green curve. 



red line of Fig. 1111 For very small values of the ratio Ao/e, 
the drapery function T) is equal to 3/4 and the difference 51 
diverges like 



51 - 



4 



e 



(59) 



Because the energy is measured with a limited accuracy, the 
relevant quantity is the averaged value of the difference 51 
which we define as 



5Ial 



1 



dE (5/(Ad) 



(60) 



This integral is performed on the energy bin that contains 
the line at Es ~ m^. Close to the line, i.e. for E ~ Es, 
the positron diffusion length Ad is given by expression (|39|) 
which we combine with relation (|59p to get 



51.. ^ M 



3.77 kpc 



^0.65 



(61) 



For a 100 GeV DM species, this leads to a difference between 
the fine-grained and coarse-grained halo integrals of 



5Iae = 5.8 X 10" 



(62) 



A 1% measurement of the energy translates into a maximal 
value of the difference (5Iab of 0.058 ~ 6%. We therefore 
conclude that even in the most extreme situation where the 
Earth lies in the region of highest DM density, we find no 
difference between the results of Section[3]and those derived 



with the fine-grained distribution (|43|) once the energy is 
properly averaged. 



5 CONCLUSIONS 

We have studied whether DM caustics in the halo of the 
Milky Way can amplify the fiux of cosmic ray antipro- 
tons and positrons which is received at the Earth. We have 
used the secondary infall model for our halo which natu- 
rally includes the caustics and assumed that the galactic 
DM is made of weakly interacting massive species. We have 
then taken into account the smearing of the caustics due 
to a present-day velocity dispersion of these particles of 
0.03 cm s~\ 

The cosmic ray antiprotons and positrons that are de- 
tected at the Earth originate from a region whose typical 
size is much larger than the shell thickness or even the shell 
separation. The associated horizon probes a large portion of 
the Milky Way and the coarse-grained average density (|30|l 
provides an adequate description. In the solar neighbour- 
hood, the coarse-grained caustic density is the same as the 
smooth NFW or isothermal cored distributions usually as- 
sumed in the literature. The difference lies at the centre of 
the Milky Way. The lower-envelope density diverges there 
with an index of ~ 2.11, hence a DM profile steeper than 
even in the NFW case. The coarse-grained shell density is 
assumed to be constant inside a sphere of radius rcut-off 
whose value is unknown. The smaller this cut-off radius, the 
more abundant DM at the galactic centre and consequently 
the stronger its signal should it reach the Earth. 

We have then computed the antiproton fiux which the 
coarse-grained shell density yields at the Earth. The reach 
of the antiproton sphere depends on the cosmic ray propa- 
gation model but is always of order a few kiloparsecs. The 
MIN configuration is associated to a rather small antiproton 
range. The associated signal is not different from what is de- 
rived assuming NFW or isothermal cored distributions. The 
MAX propagation model is characterized on the contrary 
by efficient diffusion taking over a moderate galactic con- 
vection. The antiproton sphere probes the inner and denser 
regions of the Milky Way. We find that the antiproton sig- 
nal is enhanced by a factor of ~ 30 should the conventional 
smooth NFW profile be replaced by the coarse-grained shell 
density (|30|) for which a cut-off radius of 300 pc has been 
assumed. The MED set of propagation parameters corre- 
sponds to the best fit to the B/C data and features the 
intermediate situation. It leads to the exciting possibility 
that the antiproton signal is only boosted above 10 GeV in 
the presence of caustics. Depending on the WIMP annihi- 
lation cross section, the antiproton fiux could be severely 
distorted at high energy as shown in Fig. [TO]where no artifi- 
cial boost factor is required. Thus a promising window opens 
up around a few hundreds of Gev where future antiproton 
measurements are eagerly awaited. 

We are less optimistic for the positron signal. Depend- 
ing on the cosmic ray propagation model, the positron fiux 
at the Earth may be enhanced in the presence of shells with 
respect to a smooth NFW or isothermal cored DM profile. 
However, this situation arises only at low energy where the 
observations are already well explained by the sole secondary 
background component. Thus DM caustics cannot provide 
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an explanation for the HEAT excess reported above ~ 10 
GeV and consequently produced in the solar neighbourhood. 
The solution so far invoked is based on a WIMP with a hard 
positron annihilati on spectrum like a Kaluza-Klein particle 
llCheng et al.|[200^ or a neutr a lino with a dominant W^W~ 
channel l|Hooper fc Silkll2005l : Ibelahave et al.ll2007h . Boost 
factors of ~ 10 with respect to a smooth NFW DM halo are 
nevertheless necessar y. Such a value is marginally possible 
(|Lavalle et al.|[2007bl ) in scenarios inspired by the A-CDM 
N-body numerical simulations which point towards the ex- 
istence of numerous and dense DM clumps inside which 
WIMP annihilation can be enhanced. We showed that the 
coarse-grained shell density (|30p does not provide an alter- 
native to explaining the HEAT excess. However, further data 
with higher precision and also a better understanding of the 
secondary positron background are needed to firmly exclude 
caustics as a possible explanation of the HEAT measure- 
ments. It remains very unlikely nevertheless that the present 
spectral form could be reproduced by caustics, at least by 
its coarse-grained distribution (|30p . 

We finally explored the possibility that the fine-grained 
distribution (|43|l could yield a strong positron flux should 
the Earth be embedded inside the densest part of a caus- 
tic, a region where the DM density reaches its peak value 
Pk, max- At high energy, the positron sphere shrinks. Averag- 
ing the fine-grained shell density (|13p by its coarse-grained 
approximation (|30p is no longer possible. Positrons that are 
received at the Earth have short diffusion lengths Ad and 
hence can be excellent tracers of nearby caustics. We inves- 
tigated here the case of a positron line. Positrons that are 
detected exactly at the line energy E — Es = have van- 
ishing diffusion length and if the Earth sits exactly inside the 
shell, an enhancement of the positron signal by a factor of 
~ 13, 000 is naively expected. However, energy is measured 
with a limited accuracy and the energy bin of the line has 
a non-vanishing width. Averaging Ad over the line bin leads 
to a diffusion length which is still far larger than the typical 
caustic thickness or even separation. Our study showed no 
difference between the results of Section[3]and those derived 
with the fine-grained distribution (|43p once energy is prop- 
erly averaged. We hypothesize that an extremely high energy 
resolution, presently unavailable, would allow to detect the 
nearby caustics. 

A word of caution is necessary at this stage though. The 
analysis of Section |4] is based on the assumption that cosmic 
rays diffuse on the inhomogeneities of the galactic magnetic 
field. The mean free path Afree of their random walk may be 
derived from the space diffusion coefficient K through the 
canonical relation 

K{E) = i Afree /3 . (63) 

In this diffusion scheme and with the MED set of parame- 
ters, positrons with energy e cover on average a distance 

Afroo = 1.1 X 10"* kpc e°-^ (64) 

before their next scattering on Alfven waves. Our treatment 
of cosmic ray propagation is definitely supported by the fact 
that Afreo is much smaller than the horizon size set by Ad 
and equation (|39p . However, if the nearest DM caustic is 
very close to the Earth and lies at a distance which does 
not exceed Afroc or if the Earth is embedded inside the dens- 
est part of a shell, the diffusion hypothesis breaks down. 



Positrons will essentially drift along the lines of the mag- 
netic field without encountering on their way any obstacle. 
Depending on the relative orientation of the local magnetic 
field with respect to the Earth and its nearest caustic, we 
could be possibly exposed to an intense high-energy positron 
fiux whose evaluation is clearly beyond the scope of this ar- 
ticle. 
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